Exciton-phonon coupling efficiency in CdSe quantum dots 
embedded in ZnSe nanowires 



S. Bounouar 1,2 , C. Morchutt 2 , M. Elouneg-Jamroz 1,2 , L. Besombes 1 , R. Andre 1 , E. 
Bellet-Amalric 2 , C. Bougerol 1 , M. Den Hertog 3 , K. Kheng 2 , S. Tatarenko 1 , and J. Ph. Poizat 1 
1 CEA-C'NRS-UJF group 'Nanophysique et Semiconducteurs', 
Institut Neel, CNRS - Universite Joseph Fourier, 38042 Grenoble, France, 
2 CEA-CNRS-UJF group 'Nanophysique et Semiconducteurs', CEA/INAC/SP2M, 38054 Grenoble, France, 
3 Institut Neel, CNRS - Universite Joseph Fourier, 38042 Grenoble, France, 

Exciton luminescence of a CdSe quantum dot (QD) inserted in a ZnSe nanowire is strongly 
influenced by the dark exciton states. Because of the small size of these QDs (2-5nm), exchange 
interaction between hole and electron is highly enhanced and we measured large energy splitting 
between bright and dark exciton states (AE € [4, 9.2] meV) and large spin flip rates between these 
states. Statistics on many QDs showed that this splitting depends on the QD size. Moreover, 
we measured an increase of the spin flip rate to the dark states with increasing energy splitting. 
We explain this observation with a model taking into account the fact that the exciton-phonon 
interaction depends on the bright to dark exciton energy splitting as well as on the size and shape of 
the exciton wave function. It also has consequences on the exciton line intensity at high temperature. 

PACS numbers: 78.67.Lt, 78.55.Et 



I. INTRODUCTION 

Semiconductor nanowires (NWs) have attracted great 
attention in the last few years since they hold great 
promise to become building blocks in tomorrow's 
nanoscale devices and circuits with vastpotential appli- 
cations ranging from nanoelectronics [l|-[3| , optoelectron- 
ics (light emitting diodes @,JB|, nanolasers [||), thermo- 
electrical energy conversion [7|, to biological or chemical 
sensors Q . Most of the N W growth methods allow for the 
variation of the chemical composition [9|, [lOj or doping 
PH along the longitudinal or radial directions. This en- 
ables the fabrication of semiconductor heterostructures, 
and more specifically of quantum dots (QDs). 

Single QDs have turned out to be excellent candi- 
dates for stable and efficient single photon sources 
[l6| . Within this category, QDs embedded in NWs have 
already demonstrated single photon emission [17j, even 
at relatively high temperature (T = 220K) [18|]. Quest 
for efficient and eventually room temperature QD sin- 
gle photon sources requires a good understanding of the 
excitonic dynamics in such systems. 

In a QD, exciton states are split by the electron hole 
exchange interaction into higher energy bright exciton 
states (BS) and lower energy dark states (DS) pj| with 
respective angular momentum of J z = ±1 and J z = ±2. 
In the present system, the energy splitting is relatively 
large, and has already been measured around AE = 6 
meV [2(| • As exchange interaction in semi-conductor ma- 
terials is proportional to the spatial overlap between the 
electron and hole wave functions, it is strongly enhanced 
in low dimensional objects. Values of AE £ [2,4]meV 
have been measured in quantum wells |2ll | and it has 
been demonstrated that they can be even higher in quan- 
tum wires [22] . Extremely large splitting were calculated 
and observed [H] in colloidal QDs (AE € [2,20] meV). 



This is the result of very good confinement of the car- 
riers owing to the small size of the QDs. Effects of the 
dark states on the QD luminescence become noticeable 
when AE 3> &bT, with hp the Boltzmann constant and 
T the temperature (20l. [241 [25| . We show that under non 
resonant pumping, the BS excitonic population leaks to- 
wards the DS that recombines non-radiatively, leading to 
a reduced excitonic light emission compared to the biex- 
citonic one. Transitions between the BS and DS states 
are due to hole spin flips assisted by phonons [26]. The 
speed of this processes depends essentially on the effi- 
ciency of the hole interaction with the phonon reservoir 
experienced by the QD (27| . 

In this paper, we show that the QD size can highly 
influence the QD phonon spectral density and modify the 
spin flip rates between the bright and dark exciton states. 
It not only has consequences on the exciton intensity at 
T = 4K but also at high temperatures. 

The paper is organized as follows. The sample prepara- 
tion and the setup are presented in section [II] In section 
imi the BS to DS energy splitting AE of several QDs 
is extracted using temperature dependent lifetime mea- 
surements. In section IIVI we present the experimental 
dependence of the spin flip rate versus AE, and suggest 
an explanation based on a theoretical model. These re- 
sults are then used in section [V] to discuss the exciton to 
biexciton line intensity ratio as a function of AE and of 
the temperature. 

II. SAMPLE AND SETUP 

The nanowires are grown by molecular beam epitaxy 
using the vapor- liquid-solid technique [28] . The substrate 
is GaAs 100 with a ZnSe buffer. After dewetting of a 
thin layer of gold at 500°C, the growth is performed at 
410°C. The nanowires have a diameter of 10 nm and a 
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length of 400 nm. The wire diameter (around 10 nm) is 
of the order of the bulk exciton Bohr diameter in CdSe 
(llnm), which means that the carriers in the QD are in 
the strong confinement regime. In order to study single 
QD, the nanowires are broken and dispersed on a sili- 
con substrate by direct contact. High resolution trans- 
mission electron microscope (TEM) images revealed very 
small QD sizes, from 2 to 5 nm in the nanowire direction 
[29l | . Other TEM experiments coupled to spectroscopy 
on other samples showed that this growth technique of- 
ten led to ZnSe encapsulated QDs. Photoluminescence 
on as-grown samples is centered around 2,35eV with a 
dispersion of ±0.07eV. 

The experimental apparatus is a time resolved mi- 
crophotoluminescence experiment setup. The excitation 
source is a frequency doubled picosecond Ti sapphire 
laser emitting at a wavelength of A = 950/2 = 475 nm 
(ie 2.6 eV). Nanowires luminescence is detected through a 
SX = 0.04 nm (ie SE = 0.2 meV) resolution spectrometer 
on a charged coupled device (CCD) camera for spectra 
or on a low jitter (40ps) avalanche photodiode (APD) for 
time resolved measurements (60 ps time resolution for 
the whole set-up). 



III. QD SIZE EFFECT ON THE EXCHANGE 
INTERACTION 
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Figure QJa) shows typical spectrum obtained on a sin- 
gle nanowire. The two transitions correspond to exciton 
and biexciton (respectively denoted X and XX ) an d were 
identified using cross-correlation techniques [20]. The 
biexciton binding energy is measured around 20 meV. 
The four relevant states of a neutral QD are represented 
on figure[IJc). The exciton is splitted in two states (noted 
DS and BS) linked by spin flip transition rates noted 7 sp i 
and 7sp2- The quantities Tx and Txx are the radia- 
tive recombination rates of the bright exciton and the 
biexciton, and Tnr is the non radiative decay rate of 
the dark exciton. These quantities have been measured 
at T = 4 K on nine QDs as 1/T X = 0.50 ± 0.05 ns, 
l/Txx = 0.30 ±0.05 ns, and 1/T nr = 5.0 ±0.5 ns, where 
the error is the dispersion amongst QDs. An example 
of power dependence under pulsed excitation is shown in 
figure Hlb) evidencing the linear and quadratic behavior 
of the X and XX line intensity respectively. In all the in- 
vestigated QDs, saturating intensity of the exciton emis- 
sion is smaller than the biexciton, which is a signature 
of the strong influence of the dark exciton [24j. When 
AE 3> fcsT, dark excitons cannot transit back to the 
bright state as there are no phonons available for such a 
spin flip and dark excitons are stored until recombining 
non radiatively. As a result, the photoluminescence of 
the exciton is less intense than that of the biexciton at 
saturation. The biexciton luminescence is not affected 
by the presence of the DS, whereas the exciton BS has 
a large probability to decay to the DS and not produce 
any photon. As spin flip rates are here of same order of 



Figure 1: (Color online) (a) Spectrum of a neutral QD at 
T = 4 K. (b) Power dependence of the X and XX lines under 
pulsed excitation at T — 4 K. (c) Level scheme and transition 
rates of a neutral QD. 
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Figure 2: (Color online) Transition rates 7 sp i and -y 3 p2 versus 
temperature. The dashed lines are the plot from the the- 
oretical model following phonon population evolution with 
temperature. 



magnitude or larger than the BS radiative decay, exciton 
intensity is generally small compared to the biexciton. 

Temperature dependence of exciton decay time give 
access to the energy splitting AE between DS and BS 
(20L l25l | . Spin flip transition rates depend on the spec- 
tral density of the acoustic phonons at the energy of the 
transition AE and on the phonon population that fol- 
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Figure 3: (Color online) Increase of splitting energies AE 
between DS and BS measured on nine QDs as a function of 
their excitonic emission energy. 



lows a Bose- Einstein distribution Nb{AE,T) = 1/(1 — 
exp(A-B/fesT)). The spin flip rates 7 sp i (bright to dark) 
and 7 sp 2 (dark to bright) can be written as: 

Tapi cx (7V B (AE,T) + l)i?(A£), (1) 
7sp2 cx N B (AE,T)R(AE). (2) 

The quantity R(AE) depends on the BS-DS energy split- 
ting AE but also on the spatial extension and shape of 
the wave function of the exciton|3(X |3l|. 

At T = OK, Nb = 0, so that the spectral density 
R(AE) cx -f sp i(0K) represents the bright to dark state 
transition rate and "f S p2 = 0. When the temperature is 
increased so that fc#T ~ AE, the bright exciton state is 
repopulated. For each temperature the transition rates 
7s P i and 7sp 2 are extracted from the time resolved lu- 
minescence and their evolution is fitted using the model 
based on the level scheme of figure Q](c)) with AE as the 
fitting parameter. 

This process is repeated on nine QDs. As shown on 
figure [3J measured splittings range from 4.2 meV to 9.2 
meV. A clear trend is appearing: splittings measured for 
low energy emission QDs are generally smaller compared 
to splittings measured for high energy emission QDs. The 
latters correspond to small size QDs, where electron and 
hole wave functions overlap very well, so that their ex- 
change energy is large, leading to large AE. Such re- 
lationships between QD sizes and energy splittings have 
been calculated for colloidal QDs [32] . Change in compo- 
sition of the CdzZni-zSe QD can also affect the emission 
energy and the value of AE. Our measurements do not 
follow exactly a smooth law. The shape of the confine- 
ment has also a strong influence on the wave function 
forms, and consequently on their correlation function. 
For example the prolate or oblate nature of the QD geom- 
etry appears to have an important effect. This explains 
why we obtained a cloud of experimental points following 
a general trend instead of a strict dependence. 



Figure 4: (Color online) Transition rate from BS to DS 7 sp i 
at 4K plotted as a function of splitting energy AE between 
DS and BS for nine differents QDs. The dotted line is a guide 
for the eyes. 



IV. EFFICIENCY OF THE EXCITON-PHONON 
COUPLING 

In figure |4j we have plotted the transition rates from 
BS to DS measured at 4K, 7 spl (4i4r), versus energy split- 
tings AE for nine different QDs. A clear enhancement 
of the transition rate for large AE can be observed. It 
can be also noted that the measured rates are comparable 
and even large than the excitonic radiative rate [24] . This 
is why the bright exciton is so depopulated compared to 
the biexciton in all the QDs investigated. 

The results displayed in figure 2] indicate that large 
AE splittings lead to extremely fast depopulation of the 
bright state suggesting that exciton-phonon interaction 
becomes more efficient with large AE. As explained, 
larger splittings correspond to small QDs. Exciton cou- 
pling to phonons not only depends on the energy needed 
for the spin flip but also on the QD size. But size and 
AE are not independent parameters and induce oppo- 
site effects on the transition rates. In order to explain 
this non-trivial behavior, it is necessary to calculate the 
phonon spectral density as a function of QD size and 
energy splitting AE for each particular QD to obtain a 
general trend. 

We shall only consider longitudinal acoustical (LA) 
compression modes since longitudinal optical (LO) 
phonons have energies of about 30 meV far above the 
dark - bright exciton splitting AE. In the following we 
will perform this calculation both in the 3D and in the 
ID case for the available phonon modes. The NW geome- 
try is obviously somewhere in-between these two extreme 
cases that will only give us a qualitative behavior. 

We first consider the 3D situation. The phonons dis- 
persion is bulk-like, approximated by the Debye law 
w(k) = cik, with ci the sound speed for LA phonons 
in the semiconductor material. Piezoelectric interactions 
and Frohlich longitudinal optical (LO) phonon couplings 
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are neglected. The exciton-phonon interaction is dom- 
inated by the hole-phonon interaction |26j. We there- 
fore consider only the latter, whose Hamiltonian can be 
written, in the second quantization representation with 
respect to the carrier states, as: 



Hh = -7= ^2 (4i a nfh,nn'(k)(bk + &-fe)> 



(3) 



knn' 



where (v n and a n are the hole creation and annihilation 
operators, bt and bk are LA phonon creation and annihi- 
lation operators. The index n represents the excitation 
level of the hole, and k the phonon mode. The coupling 
constant is defined as 



fh,nn'(k) = ah 



1 hk 

2pVci 



F nn ,{k), 



(4) 



where ah is the deformation potential for holes, V is the 
unit cell volume, p is the CdSe volumic mass. The quan- 
tity F nn > (k) is a purely geometrical form factor given by 



F nn ,(k) 



d J r^(r)e Ifc >„-(r) 



(5) 



In order to evaluate the coupling constant for the low- 
est hole state fn{E) we consider the harmonic oscillator 
potential ground state as the wave function of the hole: 



ip(r) 



1 



exp 



(G) 



where r± is the position component in the xy plane and 
I _l , lz are respectively the in plane, and out of plane (z 
direction) localization widths. For this wave function the 
form factor is easily found as 



F n (k) = exp 



kilj 



k z l z 



(7) 



For the lowest hole state, the phonon spectral density in 
the QD is: 

R(E)=±(N B (E) + 1) 

x jf E F " (*) F n ( fc ) W E - E ^ + 5 ( E + E ^- 

k 

(8) 

After performing the summation over k in the contin- 
uum limit, introducing the quadratic density of state of 
the phonons corresponding to the 3D case, the phonon 
spectral density is : 



with 



R(E) = R Q E 3 g(E), 

K) 2 



Ro 



8ir 2 hpcf 



(9) 



(10) 



The quantity R a contains all material parameters. The 
cubic dependence is due to the quadratic phonon density 
of state, and the function g(E) is a function of the energy 
and of the geometrical parameters of the QD [33] : 



1(E) 



C cos ( exp 



2h 2 c 2 



x cos % + 



: sin C 



(11) 



where £ is the angle of the wave vector k with respect to 
the normal to the z direction. 

Describing coupling of the hole to phonons confined 
in a nanowire the same way as in a 3D semi conductor 
bulk matrix seems a rough approximation. So we also 
considered the ID case in which the nanowire is taken 
as an infinitely thin monomode wave guide. The phonon 
density of state is constant and we consider that only the 
phonons propagating along the nanowire (z direction) can 
couple to the hole whose wavefunction is taken as iph (r) cx 
exp[— (l/2)(z/l z ) 2 ]. The spectral density becomes: 



R lD (E) cx Eg 1D (E), 



with 



9id{E) = exp 



HE 2 



2h 2 c 



(12) 



(13) 



The phonon spectral density is linear with energy but the 
geometrical factor gio(E) has the same gaussian energy 
dependence as g(E). 

Since at T = 4K, 7 sp i(4i^) cx R(AE), the BS to 
DS spin flip rate evolution from dot to dot can be de- 
scribed by calculating the phonon spectral density for 
each QD. Three calculated 3D phonon spectral densities 
corresponding to three different l z values (l z = 2,4,6 
nm) are plotted in figure (5][a) . Because of its increase 
with energy (linear for ID, cubic for 3D), R(AE) cor- 
responding to QDs with smaller l z have their maximum 
enhanced and shifted toward the higher energies. As a 
consequence, in small QDs, high energy phonons couple 
more efficiently with the hole. We propose to explain 
qualitatively the increase of the transition rate observed 
in figure @] by this size effect. 

The hole-phonon coupling efficiency, (AE) depends 
on the QD size l z and on the energy splitting AE. We 
determine the relation between l z and AE by evaluating 
the short range hole-electron exchange energy in the QD: 



AE = AE 



■■id 



1 



|^ 3rf (0)|^ 



\Mr)\ Z \Mr)\ Z d\ (14) 



where ip e /h(r) are the electron/hole wave functions, 
AE 3d is the exchange energy in the bulk material ZnCdSe 
[H (AE 3d = 0.19 meV, for a Cdo.5Z1io.5Se composition 
of the QD as measured in high resolution TEM exper- 
iments jl!), and |<^ 3ti (0)| 2 = with a* B the Bohr 
radius of the free exciton. To match the measured val- 
ues of AE (from 4 meV to 9 meV), we set the lateral 
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Figure 5: (Color online) (a) 3D Phonon spectral density for 
3 different QD sizes along the z direction (red curve: l z =2 
nm, green curve: Z z =4 nm, blue curve: l z =6 nm). Smaller 
QD z dimension shifts the maximum of the coupling constant 
toward the higher energies, and enhances its relative value, 
(b) Calculated 3D hole-phonon coupling efficiency vs tran- 
sition energy AE. (c) Calculated ID hole-phonon coupling 
efficiency vs transition energy AE. 
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Figure 6: (Color online) Calculated saturation ratio Ix/Ixx 
under pulsed excitation as a function of AE and temperature. 
The color legend represents Ix/Ixx saturation ratio. The cold 
colors show low exciton intensity and the hot color high in- 
tensity. 



confinement parameter l± at 4nm. With l z linked to a 
corresponding AE we make l z vary from 2 nm to 6 nm 
and we calculate Ri z {AE). 

The results, for the 3D and ID cases described above, 
are plotted in figure[5](b) and (c). Because of the simplic- 
ity of the considered exciton wave function, the aim of the 
calculation is neither to fit the experimental data nor to 
obtain a quantitative estimation of the transition rates. 
However, we can show that the hole-phonon coupling ef- 
ficiency is increasing with AE in both cases, which is not 
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Figure 7: (Color online) Experimental exciton-biexciton sat- 
uration ratio under pulsed excitation versus temperature for 
two QDs with respectively AE = 4.5 meV (blue squares) and 
AE = 8.5 meV (red circles). Dashed lines are vertical cut 
(fixed AE) from the model used in figure [6] 



a trivial result as the cut-off imposed by the dimensions 
of the QD makes this efficiency vanish for higher energy 
phonons. The effect of the confinement dimensions on 
the phonon spectral density can be a good explanation 
for the trend measured by the experiment in figure[H The 
two situations considered here are extreme cases and we 
can expect a real nanowire geometry to impose an inter- 
mediate behavior for the hole phonon coupling. 



V. EXCITON LUMINESCENCE INTENSITY 

This increase of the exciton-phonon coupling efficiency 
with energy splitting has some consequences on the tem- 
perature dependent exciton-biexciton saturation ratio 
under pulsed excitation. This calculated ratio is repre- 
sented in figure[5]as a function of temperature and energy 
splitting AE. We used the level scheme presented in fig- 
ure [T] and the values of the transition rates given in sec- 
tion lllll The relation between ^^(AK) and AE is taken 
from the function used as a guide for the eyes in figure [5] 
We also considered the increase of the dark exciton non 
radiative recombination r nr (T) with temperature follow- 
ing an Arrhenius law T nr (T) ~ exp(—E a /kBT), with an 
activation energy E a . For the two QDs studied in fig. [7] 
we have measured E a = 30 ± 5 meV and we have taken 
E a = 30 meV in the model used in figures [5] and [7J The 
particularity here is that the DS lifetime at T — 4 K 
(1/T nr = 5 ns) is very short compared to values reported 
in other systems (up to lfis [25]) and is of the same or- 
der of magnitude as the exciton radiative lifetime when 
temperature is raised up to only several tens of K. A 
characteristic temperature behavior is shown in figure [5] 
The exciton intensity is very small at low temperature 
(particularly in the high splitting region), increases as 
temperature is increased up to T = 50 K, and finally 
decreases with higher temperatures. The effect of repop- 
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ulation of the bright exciton due to the DS to BS spin flip 
is compensated and overwhelmed by the exciton popula- 
tion loss caused by the faster dark exciton non radiative 
recombination. As shown in figure representing the 
exciton-biexciton ratio versus temperature for 2 different 
QDs with different measured AE (4.5 meV and 8.5 meV 
) , this exciton line intensity decrease is more sensitive for 
large AE. However, at 300K, the bright exciton is still 
less luminescent for large splittings. It can be noted that 
the model fits rather well the experimental data. 



VI. CONCLUSION 

In summary, the exciton-phonon coupling efficiency is 
highly influenced by the QD size. This can explain the 
observed enhancement of the bright to dark spin flip rate 
with increasing splitting energy. As a result, the satura- 
tion intensity of exciton transition is a lot weaker than 
the biexciton one. This effect is all the more important 
as the QD exchange splitting AE is larger. The dom- 
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